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ABSTRACT 


Complex  constitutive  system  parameters  are  used  to  derive  the  form  of  Lagrange's 
equations  for  a  non-conservative  continuous  system  with  a  steady-state  sinusoidal  excitation. 
The  resulting  form  of  Lagrange's  equations,  for  the  case  of  a  dissipative  uniform  transmission 
line,  is  shown  to  yield  the  well-known  results  derivable  from  the  telegrapher's  equations.  The 
concept  of  time-averaged  Lagrangian  density  is  introduced,  including  the  correct  form  for  the 
time-averaged  energy  density  dissipation  function. 
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PREFACE 


High-frequency  (HF)  antenna  elements  for  radar  and  communications  systems  often  consist 
of  some  form  of  vertical  monopole  element  with  a  ground  plane  resting  on  or  in  close  proximity  to 
the  earth.  The  ground  plane  is  usually  made  as  large  as  is  economically  feasible  to  alleviate,  at 
large  ranges,  the  reduction  (caused  by  the  earth's  multipath  pattern  null  on  the  horizon)  of  low- 
angle  directive  gain,  and  the  reduction  (caused  by  power  absorption  and  dielectric  discontinuity  at 
the  earth  interface)  of  the  fraction  of  available  transmitter  power  that  is  radiated  into  the 
hemisphere  above  the  ground. 

The  objective  of  MITRE  sponsored  research  project  91260  is  to  develop  computer 
programs  and  accurate  models  for  predicting  the  element  pattern,  radiation  efficiency,  and  input 
impedance  of  HF  monopole  elements  on  ground  planes  of  various  forms  that  rest  on  or  are  in 
close  proximity  to  earth.  The  results  will  provide  better  analytical  capability  to  optimize  the 
performance  and  minimize  the  cost  of  the  antenna  elements  and  their  associated  ground  planes 
for  HF  antenna  arrays.  The  results  could  also  support  a  subsequent  experimental  program 
investigating  the  peformance  of  monopole  elements  under  nonideal  operating  and 
environmental  conditions. 

A  ground  plane  of  interest  is  a  wire-mesh  screen  of  arbitrary  grid  geometry.  One  approach 
to  determining  the  steady-state  wave  propagation  along  a  mesh  screen  above  earth  is  a  Floquet 
mode  analysis  as  illustrated  by  Hill  and  Wait  for  a  rectangular  bonded  wire  mesh  geometry  [1], 
Another  approach  is  to  consider  the  ground  screen  and  earth  as  a  continuous,  dissipative  non- 
uniform  transmission  line  with  the  mesh  screen  as  one  conductor  and  the  earth  as  the  return 
conductor.  The  line's  series  inductance  per  unit  length,  series  resistivity,  shunt  capacitance  per 
unit  length,  and  shunt  conductivity,  as  a  function  of  radial  position  along  the  screen,  are  usually 
determinable  from  the  mesh  geometry  and  complex  permittivity  of  the  earth. 
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For  the  latter  approach  of  a  dissipative  non-uniform  transmission  line,  the  propagation 
constant  may  be  determined  either  from  Maxwell's  equations  (together  with  Kirchhoff's  voltage 
and  current  laws)  or  from  Lagrange’s  equations.  Maxwell’s  equations  are  convenient  if  a  suitable 
coordinate  system  can  be  found  to  express  the  fields.  Lagrange's  equations  are  convenient  if  the 
stored  energy  densities  and  energy  dissipation  function  are  readily  expressible  independently  of 
the  coordinate  system. 

For  example,  it  has  been  suggested  that  Lagrange's  equations  might  be  advantageous  for 
a  distributed  transmission  line  in  which  there  is  mutual  coupling  among  elements  such  as  in  an 
infinite  array  of  closely-packed  monopole  elements  on  a  mesh  groundscreen  above  earth. 
However,  before  Lagrange's  equations  can  be  applied  to  such  a  case,  one  must  determine  the 
appropriate  form  of  Lagrange's  equations.  Although  the  forms  of  Lagrange's  equations  for 
dissipative  discrete  systems  and  non-dissipative  continuous  systems  are  well  known,  the  form  of 
Lagrange’s  equations  for  dissipative  continuous  systems  are  not  readily  found  in  the  literature. 
Furthermore,  most  of  the  applications  of  Lagrange’s  equations  are  for  transient  excitations 
whereas  the  form  of  Lagrange’s  equations  for  steady-state  excitations  is  not  as  readily  found  in 
the  literature.  It  is  also  not  clear  whether  the  use  of  Lagrange’s  equations  for  such  an  application 
will  offer  any  particular  advantage  over  approaches  based  on  Maxwell's  equations. 

This  paper  attempts  to  address  these  issues  by  applying  Lagrange's  equations  to  the  case 
of  a  continuous,  dissipative  uniform  transmission  line  with  a  steady-state  excitation. 
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SECTION  1 


INTRODUCTION 


Continuously-distributed  dissipative  uniform  transmission  lines  have  been  employed 
in  communication  and  radar  systems  since  the  earliest  conception  of  such  systems. 
Examples  of  such  lines  are  the  parallel  two-wire  line  and  the  coaxial  line.  Such  lines  are 
usually  analyzed,  for  sinusoidal  steady-state  excitation  and  for  the  condition  that  the 
radiation  from  the  lines  is  negligible,  by  use  of  the  telegrapher's  equations. 

The  purpose  of  this  paper  is  to  present  a  form  of  Lagrange's  equations  and  the 
associated  Lagrangian  density,  for  a  dissipative  uniform  transmission  line  with  steady-state 
excitation,  that  yield  identical  results  as  the  telegrapher's  equations.  One  motivation  for 
doing  so  is  to  provide  a  foundation  for  future  effort  in  applying  Lagrange’s  equations  to 
more  complex  transmission  lines,  as  discussed  in  the  Preface.  A  second  motivation  is  to 
evaluate  whether  Lagrange's  equations  offer  any  particular  advantage  for  such  an 
application. 

The  forms  of  Lagrange's  equations,  for  transient  excitation  of  dissipative  discrete 
systems  and  non-dissipative  continuous  systems  are  well  known.  However  the  form  of 
Lagrange's  equations  for  sinusoidal  steadv-state  excitation  of  dissipative  continuous 
systems  is  not  readily  found  in  the  literature.  This  paper  attempts  to  address  this 
deficiency. 

Some  of  the  concepts  associated  with  Lagrange's  equations  that  are  addressed  in  this 
paper  are: 

( 1 )  the  form  of  Lagrange’s  equations  for  distributed  vs.  discrete  systems; 


(2)  the  time-averaged  Lagrangian  density  for  a  continuous  system  with  a  force 
field  derivable  from  a  generalized  potential  comprising  scalar  and  vector 
potentials; 

(3)  the  time-averaged  energy  density  dissipation  function  for  a  non-conservative 
electromagnetic  system  with  series  and  shunt  ohmic  losses; 

(4)  the  relationship  of  Lagrange's  equations  to  Hamilton's  principle  and  the 
calculus  of  variations;  and 

(5)  the  definition  in  the  calculus  of  variations  of  the  partial  derivative  of  a  non- 
analytic  complex  function. 

The  well  known  results,  for  the  propagation  constant  of  a  continuously-distributed 
dissipative  uniform  transmission  line  with  sinusoidal  steady-state  excitation,  are  derived 
from  the  telegrapher's  equations  in  section  2  and  from  Lagrange's  equations  in  section  3 
supported  by  an  appendix.  The  conclusions  are  given  in  section  4. 
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SECTION  2 


TELEGRAPHER'S  EQUATIONS 

Consider  an  increment  dx  of  a  dissipative  uniform  transmission  line  (see  figure  1) 
with  distributed  series  inductance  L  dx  (henries),  shunt  capacitance  C  dx  (farads) 
series  resistance  R  dx  (ohms)  and  shunt  conductance  G  dx  (Siemens)  along  the  x  axis. 

For  a  steady-state  sinusoidal  wave  of  radian  frequency  co  (rad/s),  the  voltage  v(x,t)  at 
point  x  and  time  t  may  be  expressed  in  the  form 

v(x,r)  =  Re[V(x)e;<i"]  (2-1) 

/(x,r)  =  Re[/(x)0  (2_2) 

It  is  assumed  that  the  power  radiated  by  the  transmission  line  into  the  surrounding 
free-space  medium  is  negligible. 

The  incremental  voltage  dv  and  incremental  current  di  are  given  by  [2]. 

dv(\,t)  =  v(x  +  dx,r)  -  v(x,r) 

=  -R  dx  i(x,t)-L  dx 

dt 

=  Re[-(/?  +  jcoL)I(\)  dx  eJQ“] 

=  Re[dV(x)e^']  (2.3) 


CHARGE  q  =  |Pdx 

.  3q  (x.  I) 


Figure  1 .  Increment  of  Dissipative  Uniform  Transmission  Line 


di(x,t)  -  i(x  +  dx,r)  -  i(x,t) 


=  G  dx  v(x,t)-C  dx—-’f) 

dt 

=  Re[-(G  +  jcoC)V(x)  dx  ejm‘] 

=  to[<ff(x)e*"]  (2.4) 

Equations  (2-3)  and  (2-4),  after  cancelling  the  common  real  part  designation  ana  the 
common  factor  eJbn,  yield  the  well-known  telegrapher's  equations: 

— =  ~(R  +  j(oL)I  (x) 

dx  (2-5) 

^-  =  -(G  +  jcoC)V(x) 

dx  (  2  -  6  ) 


Differentiating  equation  (2-5)  with  respect  to  x  and  then  substituting  dl/dx  from 
equation  (2-6), 

d2V(x) 

V —  -  ( P  +  jcoL)(G  +  jcoC)V(x)  =  0 
dx1  (2-7) 

Noting  that  (R  +  jouL)  (G  +  jcoC)  is  a  complex  constant,  equation  (2-7)  has  the  general 
solution 

Vfx)  =  exp  (-yx)  +  V_  exp  (yx)  (2-B) 

where  V+  and  V.  are  arbitrary  complex  constants  and 

y2  =  (R  +  jcaL)(G  +  jcvC)  (2-9) 
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Differentiating  equation  (2-8)  with  respect  to  x, 


dVjx) 

dx 


= -)V+ exp(-yx)  +  )V_exp(yx) 


Substituting  equation  (2-10)  into  equation  (2-5), 


/(x)  =  - 


dV(\)  dx 
R  +  jcoL 


y 

R  +  jcoL 


[Kexp(-yt)-V/.exp(yx)J 


=  (1/Z.)  [V+exp(-yt)-V  exp(yx)] 


(2-10) 


(2-11) 


where 

Z0=[(R  +  jcoL)/(G  +  jcoC)]U2=  characteristic  impedance  of  the  transmission 

line. 


The  quantity  Y  in  equations  (2-8)  through  (2-11)  is  defined  as  the  complex  wave 
propagation  constant  (or  simply  "propagation  constant").  Equation  (2-9)  may  be 
rewritten  as 


y2-(R  +  jcoL)(G  +  jcoC)  =  0  (2-12 

which  may  be  thought  of  as  the  characteristic  equation  for  the  propagation  constant  y. 
Equation  (2-12)  is  a  quadratic  equation  with  two  roots  for  y.  Since  the  transmission 
line  of  figure  1  is  a  passive  medium,  the  propagation  constant  y  is  that  root  of 
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equation  (2-12)  for  which  Rey  ^  0  and  ImY  ^0  for  o)>0.  The  propagation  constant 
y  =  a+  jP  is  therefore  given  by 

y  =  [(R  +  jcaL)(G  +  jcoC)]m;  a  >  0,  P>0  for  co>0  (2-1 

where 

a  =  Re  y  =  attenuation  constant 
P  =  Imy  =phase  constant  ( m _I) 
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SECTION  3 


LAGRANGE'S  EQUATIONS 

3.1  Currents,  Charges,  Charge  Densities 

For  the  sinusoidal  steady-state  the  total  charge  q(x,t)  (coulombs)  in  the  incremental  circuit 
of  figure  1  may  be  expressed  in  the  form 

<7(x,r)  =  Re[<2(x)e;a"]  =  Re[Z(x,r)]  =  (1  /  2)[Z(x.r)  +  Z(x,r)]  ^ 

where 

Z(x,r)  =  Q(x)ejC0‘  =  complex  total  charge 

Z(x,r)  =  Q(x)e~10)‘  -  complex  conjugate  total  charge 

and  the  superscript  —  denotes  complex  conjugate.  The  charge  density  p(x,t)  (coulombs/m)  is 
related  to  the  charge  q(x,t)  by 

q  =  \p  dx 

or  (3-2) 

p  =  dq  /  dx  =  qx 

The  current  i(x,t)  (amperes)  is  related  to  the  charge  q(x,t)  by 

i  =  dq/dt  =  q  (3-31 


3-  1 


where  the  superscript  •  denotes  3  /  dt .  Differentiating  equation  (3-20)  with  respect  to  t, 

p  =  qx  =  ix  =  di/ dx  (3-4) 

Equation  (3-4)  is  the  continuity  equation  for  charge  and  current.  Equation  (3-4)  differs  by  a  minus 
sign  from  the  conventional  form  of  the  continuity  equation  because  in  figure  1  the  positive  charge 
at  node  "a"  of  the  shunt  circuit  is  identified  with  positive  current  flowing  inis  the  node.  In  the 
conventional  derivation  of  the  continuity  equation,  positive  charge  is  identified  with  positive 
current  flowing  out  of  the  node  [3]  (corresponding  to  positive  current  density  at  a  surface  in  the 
direction  of  the  outward  normal  to  the  surface). 

The  differential  current  di,  in  the  shunt  circuit  consisting  of  the  capacitance/length  C  in 
parallel  with  the  conductance/length  G,  is  found  from  equations  (3-2)  through  (3-4)  to  be 


di  =  dq  =  p  dx  =  qxdx  =  ixd  x 


(3-5) 


3.2  Current  and  Voltage  Constraints 

Conservation  of  charge  and  current  requires  that  the  sum  of  the  charges  or  currents  into  a 
node  must  equal  the  sum  of  the  charges  or  currents  out  of  the  node.  Accordingly,  the  charge 
density  p,  charge  dq,  and  current  di  into  node  "a”  are  related  to  the  charge  densities,  charges, 
and  currents  in  arms  1  and  2  of  the  shunt  circuit  by 


P  =  P,  +  P2 

(3-6) 

di  =  dil  +  di2 

(3-7) 

dq  =  dqx  +  dq2 

(3-8) 

From  equation  (3-6)  it  follows  that 
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(3-9) 


P  =  P,+P2 


The  voltage  across  the  unit  capacitance  C  must  equal  the  voltage  across  the  unit 
conductance  G.  Accordingly,  from  equation  (3-5)  it  follows  that 

v(x  +  dx  r)  ~  =  p^x  =  A  =  =  P2 

Cdx  Cdx  C  C  Gdx  G  (3-10) 


Therefore, 

p2(x,t)  =  (G  /  C)pl(x,t)  (3--M) 

Substituting  equation  (3-10)  into  equation  (3-9), 

p  =  pl  +  (G/C)Pi  =  qu  +  (G/C)glx  (3-12) 

Multiplying  both  sides  of  equation  (3-12)  by  dx  and  taking  note  of  equation  (3-5),  equation  (3-12) 
reduces  to 

q  =  qx  +  {GIC)qx  (3-13) 

3.3  Time'Averaged  Lagrangian  Density  in  the  Absence  of  Ohmic  (Dissipative) 
Losses 

For  a  steady-state  sinusoidal  waveform,  the  time-averaged  Lagrangian  density  {l)  (J/m)  for 
the  circuit  of  figure  1 ,  in  the  absence  of  ohmic  losses  (R=G=0),  is  given  (with  reference  to 
equation  (A-8)  of  the  appendix)  by 
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(l(Z,Zx,Z,Z,Zx,Z, x,t)^j  =  (m)  —  (e),  R  =  G  =  0 


(3-14) 


where 

(m(Z,Z,Z,Z,x)J  =  time-averaged  stored  magnetic  (kinetic)  energy  density  (J/m) 
for  R=G=0. 

(e{ZxZ,Zx,Z,x)}  =  time-averaged  stored  electric  (potential)  energy  density  (J/m) 
for  R=G=0. 


Z(x,t)  =  complex  total  charge  defined  by  equation  (14). 

The  time-averaged  stored  magnetic  and  electric  energy  densities  for  R=G=0  are  given  by: 


(m)  = 


(^Lq2^J  =  {\ /  4)LZZ  =(1/  4)LZ,Z,,  R  =  G=Q) 


(3-15) 


(e)  =  (l  /  2)Cvz(x  +  dx,t))  =  ((1  /  2)C(qXx  /  C)2)  =  (1  /  4)(1  /  C)ZUZU,  R  =  G  =  0  (3-16) 


where 

f\<7>v  are  given  by  equations  (3-5),  (3-13),  and  (3-10),  respectively. 

Z;(x,r)  =  (?,(x)e;  _  comp|ex  charge  at  the  distributed  capacitance  C 

such  that  q,(x,0-Re[Zt(x,/)]  (3-17) 

Lagrange's  equation  for  the  circuit  of  figure  1  with  R=G=0  is  given  by  equation 
(A-10)  of  the  appendix  as: 
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(3-18) 


d_ 

dt 


f  j  / 


m 

\  3Z\  J 


d_ 

dx 


v  az.,  J  az, 


where  the  partial  derivatives  of  <£)  are  defined  in  the  appendix  from  the  calculus  of  variations  by 
treating  Z,Zt,Z,Z,Zx,Z  as  independent  variables.  The  time-averaged  Lagrangian  density  (l)  is 
found  from  equations  (3-14)  through  (3-17)  to  be 


<£)  =  (!  /4)LZ£  -  (1/  4)(1  /  C)ZlxZu,  R  =  G  =  0 


(3-19) 


It  follows  that: 


^  =  (1/4)^  =  (1/4)7^ 
oZ, 


(3-20) 


dfd(L)^ 
dt  ^  <JZX  j 


=  (1  /  4)LZ,  =  -(1  /  4)co2LZ,  =  -(1  /  4)ruzLQ, 


.2;  n  jtuc 


(3-21) 


^  =  -(1  /  4X1  /C)Z„ 

»i« 


(3-22) 


V  azu  y 


=  -(i  /  4)(i  /  c)z,„  =  -d  /  4)(i  /  c)Glwy  • 


(3-23) 


a<£) 


=  0 


(3-24) 


Substituting  equations  (3-20)  through  (3-24)  into  equation  (3-18), 
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=  0,  R  =  G  =  0 


-(1  /  4)6>2L<2,e“"  -  (1  /  4)(1  /  C)Qlxxe‘“' 

or 

(?,„  +  (o2LCQl  =  0,  R  =  G  =  0 
Designating  y2  =  -to2  LC,  equation  (3-25)  may  be  written  as: 

a„-r%  =  0>  /?  =  G  =  o 

Equation  (3-26)  has  the  general  solution: 

QM  =  &e-"  +  Q_e*\  R  =  G  =  0 
where: 

Y  =  a  +  jp  =  propagation  constant  =  jco  (LC)1/2,  R=G=0 

a=  attenuation  constant  =  0,  R=G=0 

P  =  phase  constant  =  co  (LC)1/2,  R=G=0 

Q+  Q.  =  complex  constants 

The  charge  q(x,t)  =  Re[Zi(x,t)]  is  therefore  given  by 

<7,(x,r)  =  Rc[(Q+e-rx  +  Q,er%)e,a,‘] 

=  Re[a  eja>-n-Q_e-l<0,+n] 

=  Re[&ej{a,-p')-Q_ej{-a,,+p')] 

=  Re{a^<ul''(,/V')1  - 0_e"'4+M},  R  =  G  =  0 


(3-25) 

(3-26) 

(3-27) 

(3-28) 

(3-29) 

(3-30) 


(3-31) 
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The  complex  charge  Zi(x,t)  in  a  non-dissipative  uniform  transmission  line  consists  of  two  waves 
travelling  with  a  phase  velocity  vp  =  (oV3)  =  (LC)'1^  ~  one  wave  of  constant  amplitude  Q+ 
travelling  in  the  +x  direction  and  the  other  wave  of  constant  amplitude  Q.  travelling  in  the  -x 
direction. 

3.4  Time-Averaged  Lagrangian  Density  In  the  Presence  of  Ohmic  (Dissipative) 
Losses 

The  transmission  line  of  figure  1  is  a  continuous  nonconservative  electromagnetic  system 
with  stored  conservative  electromagnetic  fields  and  dissipative  electromagnetic  forces  acting  on 
moving  charges.  The  total  force  field  (the  so-called  "Lorentz  force”)  is  derivable  from  a 
generalized  potential  that  is  a  combination  of  scalar  and  vector  magnetic  potentials  [4,5], 
Lagrange's  equation  for  this  total  force  field  may  be  written  in  the  same  form  as  equation  (3-18) 
provided  that  one  substitutes  for  the  time-averaged  Lagrangian  density  (L),  given  by  equation 
(3-19),  a  modified  time-averaged  Lagrangian  density  ( L)  (the  prime  indicates  "modified”)  that  is  a 
function  of  the  generalized  potential  [4,5]. 

The  modified  time-averaged  Lagrangian  density  ( L  )  may  be  expressed  as: 

<£.’>  =<£)+<?)  (3-32) 


where: 

(  l)  =  is  given  by  equation  (3-19) 


(  f)  =  time-averaged  energy  density  dissipation  function  (J/m) 
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The  time-averaged  energy  density  dissipation  function  ( 9)  is  generally  a  complex  quantity  [6], 
The  dissipation  function  (  7)  is  readily  determined  by  defining  a  complex  inductance/length,  L\ 
and  a  complex  capacitance/length,  C',  by: 


L'  =  L  +  (1/jco)R  =  L  -  j  (R/m) 


(3-33) 


C'  =  C  +  (1/jco)G  =  C  -  j  (G/co) 


(3-34) 


It  should  be  noted  that  L'  and  C'  have  the  dimensions  of  complex  permability  and  complex 
permittivity,  respectively,  since  [L]  =  [henries/m]  =  (permeability],  [C]  =  farads/m]  =  [permittivity], 

[R]  =  [ohms/m]  =  [resistivity],  and  [G]  =  [siemens/m]  =  [conductivity].  Substituting  L'  for  L  and  C'  for 
C  in  equation  (3-19): 

<r)  =  (m')  -  (e')  =  (1  /  4)L'Z,2,  -  (1/  4)(1  /  C)Z„Z„ 

=  a  /  4)tZ,f,  -  (1  /  4)J(K  /  <0  )iz£  -  (1  /  4)[c  -  j(G  /  ffl)]"Z„Zu  (3.35) 


where: 


<  m )  =  time-averaged  complex  magnetic  energy  density 


=  (1  /  4)L'Z1Z1  =  (1  /  4)LZ,Z,  -  (1  /  4 )j(R  /  (o)LZ,Z, 


(3-36) 


( t)  =  time-averaged  complex  potential  energy  density 


=  (1  /  4)(1  /  C)ZUZU  =  (1  /  4  )[C  -  j(G  /  co]"ZuZu 


(3-37) 


Substituting  equations  (3-35)  and  (3-19)  into  equation  (3-32),  the  time-averaged  energy  density 
dissipation  function  (  7)  is  given  by 
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(T)  =  (1  /  4  ){L'-L)ZX  ~  (1  /  4)[(1  /  C)-  (1  /  C)]Zj]z 
=  (1  /  4)(/?  /  j(0)Z.X  +  (1  /  4)[(G  /  C)  /  (yryc  +  G)]ZuZlx 


Substituting  <  L  >  for  <  L)  in  equation  (3-18),  Lagrange's  equation  for  the  system  of 
figure  1  becomes: 


f  d(L')] 

,± 

dt 

<s,  J 

dx 

dZx 


(3-39) 


where  <L)  is  given  by  equation  (3-35).  It  follows  from  equations  (3-20)  through  (3-25)  that 
equation  (3-39)  reduces  to. 


Q^xx  +  co2LCQl  =  0 


(3-40) 


which  has  the  general  solution 

=  &  exp  (-y'x)  +  Qlexp(y’x) 


(3-41) 


where: 

Y-  7 +jP'=  propagation  constant  =  jco(L'C')u 2 
=  jcd[L  +  ( R  /  jo))]u2[C  +  ( G  /  j(o)]m  =  [{jcoL  +  R){jcoC  +  G)}m 
=  jcc{LC)in{[  1  -  (JR  /  <uL)][l  -  (yG  /  otC)]}"2 

(3-42) 
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a'=  attenuation  constant  =  V2  +  (02l}){G2  +  ry2C2)]  2  +  RG  -  co2LC 

(3'=  phase  constant  =  V2j[(/?2+  o)2L2)(G2  +  cd2C2)]'/2  +  RG  +  oj2Lc| 

Q+,  Q_  =  complex  constants 


(3-43) 


(3-44) 


The  charge  qi  (x,t)  =  Re  [Zi  (x,t)]  is  therefore  given  by: 


<7i(x,0  =  Re{[Q^  exp  (-y'x)  +  Q_  exp  (y’x)]} 

=  Re  [2^.  exp  (jco  t-y'x )  -  Q_  exp  (- jco  t  +  y'  x)] 
=  Re  [Q;e-a'*e;(r'<a"/,‘)-Q.ea'V;('r<u'^',) 

=  Re  [ ’ae-a'*e>‘uHx/,,v))-aea'^';<ul'+(,/v')r 


(3-45) 


The  complex  charge  Zi(x,t)  in  a  dissipative  uniform  tranmission  line  consists  of  two  waves 
travelling  with  a  phase  velocity 

v,  =  (a>  /  P')  =  ^C)'m  {[1  -  (JR  /  ojL)][1  -  {JG  /  ojC)}\m 

-  one  wave  of  decreasing  amplitude  Q*  e~a  1  travelling  in  the  +x  direction  and  the  other  wave  of 
decreasing  amplitude  (2-  e~a  *  travelling  in  the  -x  direction. 


It  should  be  noted  that  the  propagation  constants  given  by  equations  (3-42)  and  (2-13)  are 
identical.  Therefore,  the  results  obtained  from  Lagrange's  equations  agree  with  those  obtained 
from  the  telegrapher's  equations,  as  should  be  the  case. 
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SECTION  4 


CONCLUSIONS 


Complex  constitutive  system  parameters  are  used  to  derive  the  form  of  Lagrange's 
equations  for  a  non-conservative  continuous  system  with  a  steady-state  sinusoidal  excitation. 
Confirmation  is  demonstrated  by  considering  the  case  of  a  uniform  transmission  line  with  series 
and  shunt  ohmic  losses,  sinusoidal  steady-state  excitation,  and  wave  propagation  in  a  transverse 
electromagnetic  (TEM)  mode.  The  resulting  voltages,  currents,  propagation  constant,  and 
characteristic  impedance  are  in  agreement  with  the  well-known  results  derivable  from  the 
telegrapher's  equations. 

The  forms  of  Lagrange’s  equations  for  sinusoidal  steady-state  excitation  and  for  dissipative 
continuous  systems  are  not  readily  found  in  the  literature.  (The  forms  of  Lagrange's  equations  for 
transient  excitations,  dissipative  discrete  systems,  and  non-dissipative  continuous  systems  are 
well  known.) 

In  particular,  some  of  the  concepts  that  are  addressed  in  this  paper  are: 

(1 )  The  form  of  Lagrange's  equations  for  distributed  versus  discrete  systems; 

(2)  The  time-averaged  Lagrangian  density  for  a  continuous  system  with  a  force  derivable 
from  a  generalized  potential  comprising  scalar  and  vector  potentials;  and 

(3)  The  time-averaged  energy  density  dissipative  function  for  a  non-conservative 
electromagnetic  system  with  series  and  shunt  ohmic  losses. 

For  the  case  of  a  dissipative  uniform  transmission  line,  Lagrange's  equations  offer  no 
advantages  over  that  of  the  telegrapher's  equations  or  equivalent  equations  based  on  Kirchhoff's 
voltage  and  current  laws  or  Poynting's  complex  power  theorems.  The  reason  is  that,  for  this  case, 
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at  least  as  much  effort  is  required  to  determine  the  stored  energy  densities  and  dissipation 
function  as  is  required  to  express  the  fields  directly. 

However,  having  confirmed  the  basic  soundness  of  the  Lagrangian  formulation  for  the 
analysis  of  dissipative  transmission  lines  and  having  established  the  form  of  Lagrange’s  equations 
for  the  simplest  case  of  a  uniform  transmission  line,  this  paper  is  expected  to  facilitate  the  use  of 
Lagrange's  equations  for  transmission  lines  with  more  complex  geometries,  particularly 
geometries  for  which  the  telegrapher's  equations  or  Maxwell’s  equations  may  not  be  tractable.  A 
Lagrangian  formulation  is  superior  in  cases  where  the  fields  are  not  readily  expressible  in  a 
suitable  coordinate  system  but  where  the  stored  energy  densities  and  energy  density  dissipation 
functions  are  expressible  without  reference  to  a  particular  set  of  generalized  coordinates^4). 

A  suitable  candidate  application  might  be  a  distributed  transmission  line  in  which  there  is 
mutual  coupling  among  several  elements  such  as  in  an  infinite  array  of  closely-packed  monopole 
antenna  elements  on  a  mesh  groundscreen  above  earth.  The  forms  of  Lagrange’s  equations 
developed  in  this  paper  should  be  applicable  to  such  a  transmission  line  but  the  detailed 
application  for  this  case  remains  to  be  investigated.  It  is  unclear,  however,  based  on  the  present 
analysis,  that  Lagrange's  equations  will  offer  any  particular  advantage  for  that  application. 
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APPENDIX 


LAGRANGE’S  EQUATIONS  FOR  A  CONTINUOUS  SYSTEM  WITH 
A  STEADY-STATE  SINUSOIDAL  WAVEFORM 


Consider  a  continuous  system  that  is  conservative  with  a  force  field  derivable  from  a  scalar 
potential  or  is  nonconservative  with  a  force  field  derivable  from  a  generalized  potential  comprising 
scalar  and  vector  potentials.  For  such  a  system  with  n  coupled  elements,  the  Lagrangian  density 
L  is  defined  as  [1,3]. 

L  —  m  —  e  —  L{^qvq1, — (A-1) 

where: 

m  =  kinetic  (magnetic)  energy  density  of  the  system 
e  =  potential  (electric)  energy  density  of  the  system 

qj  =  position  (charge)  of  i,h  element  that  is  a  function  of  spatial  and  time  coordinates  x 
and  t,  respectively 


superscript  •  denotes  d/dt 
subscript  x  denotes  3/3x 


Hamilton's  principle  states  that  the  motion  (field)  of  the  system  from  point 
Pi  (xi  ,ti )  to  point  ?2  (*2.t2) is  such  ,hat  the  line  inte9rat  ®  is  9iven  by  i2-3) 


O  =  f'3  f  *2  L  dx  dt  =  extremum  for  the  path  of  motion 

J/,  Jx, 
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i.e.,  the  system  will  travel  only  that  path  for  which  the  integral  of  *he  Lagrangian  density  from  point 
Pi  to  point  P2  is  either  a  minimum  or  maximum. 


Equivalently  equation  (A-2)  states  that  the  variation  80  of  the  line  integral  <j>  for  a  fixed 

W  xi.  x2  iszer0  I2’3!- '  e  - 

80  =  5 


JT 


Ld\dt  =  0 


(A-3) 


With  the  substitution  of  equation  (A-1)  into  equation  (A-3)  and  the  utilization  of  the  calculus  of 
variations  [1 ,2.3],  equation  (A-3)  reduces  to 


f’pjy 

’  d 

(is? 

d 

4_  _ 

rdL 

dL]\ 

dt 

dx 

"  ^  Jj 

dx  dt  =  0 


(A-4) 


where  the  partial  derivatives  are  defined  in  the  calculus  of  variations  by  treating  qj,  qjx,  q*  as 
independent  variables.  Since  the  qj  variables  are  independent  80  =  0  if  and  only  if  the 
coefficients  of  dqj  separately  vanish  [1,3].  Accordingly, 


didS 

dtydqiy 


'IL 

A, 


i  =  1,2, . n 


(A-5) 


The  n  equations  given  by  equation  (A-5)  are  Lagrange's  equations  for  a  continuous  system. 

For  continuous  systems  with  a  steady-state  sinusoidal  waveform,  the  time-averaged 
Lagrangian  density  may  be  more  convenient  to  compute  than  the  instantaneous  Lagrangian 
density  given  by  equation  (A-4).  Lagrange's  equations  that  are  a  function  of  a  time-averaged 
Lagrangian  density  may  be  derived  by  dividing  both  sides  of  eq  (A-3)  by  t2  -  t, .  Accordingly, 
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_J—  <5<Z>=-i— dxdt\  =  8  f-i-Vp£  dxdt  =SP(L)dx 
h-t i  t2-tx  IA  J*>  J  J  (A-6) 

where  {£)  =  — - —  J  1  £  dt  =  time  -  averaged  Lagrangian  density 


and  the  brackets  (  >  indicate  time-averaged  Lagrangian  density. 

For  a  steady-state  sinusoidal  waveform  (t2-t-|  -*»)  the  time-averaged  Lagrangian  density  is 
in  general  a  function  of  ZjZ,-,  Z^Z^,  Z,Z,  where 


=  £(x)  e;“'  with  Re[Z,(x,/)]  =  <7,(x,/)  (A_7) 

For  such  a  waveform  the  time-averaged  Lagrangian  density  <  L )  may  be  written  in  the  functional 
form 

=  y.a.Z.Z,  +  £,ZlzZu  +ciZiZj 
,=i  (A-8) 


(£)  —  ^£  Zu,  Z,,  Z„  Zu,  Z;j,  x,r 


where  q,  bj,  q  are  complex  constants.  Equations  (A-8)  and  (A-6)  are  analogous  to  equations 
(A-l)  and  (A-3),  respectively,  Lagrange's  equations  that  are  a  function  of  the  time-averaged 
Lagrangian  given  by  equation  (A-8)  are  exactly  analogous  to  equations  (A-5).  Accordingly,  the 
corresponding  Lagrange's  equations  are  given  by: 


dfdjL)}  d(d(L? 
dt{  JZ,  J  dx\dZixJ 


i  =  1,2,. ..n 


(A-9a) 
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(A-9b) 


d 

(W) 

d 

J _ 

f  d(L)} 

d(L) 

dt 

{ 4  J 

T 

dx 

dz, 

It  should  be  noted  in  equation  (A-9)  that  (  L)  is  not  an  analytic  function  of  the  complex 
parameters  Z,Z  because  <  l)  given  by  equation  (A-8)  does  not  satisfy  the  Cauchy-Riemann 
equations.  Therefore,  the  derivative  of  <  l)  with  respect  to  either  Z  or  Z  is  not  defined 
(lim  a(l)/az  as  AZ— >0  is  dependent  upon  the  path  taken  in  the  Z  plane).  However,  the  partial 
derivatives  in  equation  (A-9)  are  defined  in  the  exact  analogous  manner  as  the  partial  derivatives  in 
equation  (A-4),  i.e.,  by  treating  ZZ  ,  ZUZU,  Z-Z,  as  independent  variables. 


Equation  (A-9a)  leads  to  a  solution  for  Z,  whereas  equations  (A-9b)  leads  to  a  solution  for 
Z  Both  sets  of  equations  give  the  same  result  for  the  propagation  constant.  Since  Z,  is 

derivable  from  Z,  and  vice  versa,  equations  (A-9a)  and  (A-9b)  are  redundant.  Only  one  set  of 
equations  need  be  considered.  Selecting  equation  (A-9b),  a  non-re dundant  set  of  Lagrange's 
equations  are  given  by: 


d 

{  d(4) 

d 

(  W) 

dt 

dx 

i  =  1,2,. ..n 


(A- 1 0) 


A  more  restricted  form  of  Lagrange’s  equations  for  a  steady-state  sinusoidal  waveform 
may  be  obtained  by  considering  only  the  space-dependent  complex  amplitude  Q(x)  of  the 
variable  q(x,t)  in  equation  (A-7).  The  corresponding  Lagrangian  density  L  is  then  of  the  form 


L=  L 


\i=i 


(A-11) 


Hamilton's  principle  states  that  the  variation  of  the  line  integral 
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Ldx  for  a  fixed  x,  and  x2  is  zero,  i.e. 


=  £dx  =0 


(A- 1 2) 


The  calculus  of  variations  yields  in  a  manner  exactly  analogous  to  that  of  equation  (A-4),  a  variation 
d0  given  by 


'  d fdL" 
dx\dQlxJ 


dL 

JQ, 


i=i 


mW  =  o 


(A-13) 


where  Qw Q>Qi  are  treated  as  independent  variables.  Since  Q,  and  0,  are  independent, 

80  =  0  if.  and  only  if,  the  coefficients  of  <50,  and  <50,  separately  vanish.  Accordingly, 


d (dL\  dL 

dQuJ  XI 


{ A- 1 4) 


d_ 

dx 


(A- 1 5) 


Since  Q,  is  derivable  from  Q,  and  vice  versa,  equations  (A-14)  and  (A-15)  are  redundant. 
Therefore,  either  set  of  equations  (A-14)  or  (A-15)  is  sufficient. 
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dpL 

dx^q„ 
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